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Abstract 

The recent observations of 50 TeV gamma radiation by HEGRA have 
the potential of determining the extragalactic flux of infrared radiation. The 
fact that radiation is observed in the range between 30 and 100 TeV sets 
an upper hmit on the infrared flux, while a cutoff at E-y ~ 50 TeV hxes 
this flux with a good accuracy. If the intrinsic radiation is produced due to 
interaction of high energy protons with gas or low-energy target photons, then 
an accompaning high-energy neutrino flux is unavoidable. We calculate this 
flux and underground muon flux produced by it. The muon flux is dominated 
by muons with energies about 1 TeV and can be marginally detected by a 
1 km^ detector like an expanded AMANDA. 



The detection by HEGRA |T|] of sources of gamma radiation with energy about 
50 TeV, if confirmed, is a remarkable observation which hopefully is the first step 
in ultra high energy gamma astronomy. Besides its own importance, this radiation 
is an excellent tracer of infrared fluxes in intergalactic space. 

The existence of an upper limit of the observed energies of photons is interpreted 
in as absorption on diffuse infrared radiation (IR). This assumption looks rea- 
sonable due to the following argument. Observation by CAS A p| has not resulted 
in detection of these sources 0. This detector has much higher flux sensitivity 
than HEGRA, but also higher energy threshold {Eth = 100 TeV), which is sharp 
and makes ineffective the observations of gamma ray sources below about 50 TeV 
0]. Therefore the cutoff of the spectrum at 50 TeV must be very sharp. The known 
acceleration mechanisms and in particular shock acceleration cannot provide such 
a sharp cutoff, while absorption naturally results in an exponential cutoff. 

On the other hand, the observation of 50 TeV gamma-radiation from the HEGRA 
sources, which are typically at a distance greater than 200 Mpc, implies that the 
flux of diffuse IR radiation is lower than was theoretically estimated (see ^ and 
references therein^). 

It is interesting to note that recently the EASTOP collaboration analyzed their 
data on Mk 421, which was detected by HEGRA with 3.8a excess over background. 
At energy higher than 40 TeV EASTOP collaboration established [g the 90% c.l. 
upper limit 1.2 ■ 10^^^ cm^^s^^ close to the flux observed by HEGRA. Within the 
statistical errors there is no serious contradiction between these two measurements. 
In our analysis we concentrate on the source 0116+319 detected by HEGRA at 5.7cr 
level with flux 1.4 • 10"^^ cm'^s"^ at E > 50 TeV. 

We shall discuss some consequences of the HEGRA observations. One of them 
is the deriviation of the flux of IR radiation consistent with these observations. If 
the intrinsic flux has its origin in the interaction of high energy protons with gas 
or target photons then an accompanying high energy neutrino flux is unavoidable. 
We shall also discuss the detectability of this flux. 

Let us parametrize the density of IR photons in intergalactic space, n{e), as 

n(e) = '^(^r, (1) 



^0 \£^0' 

where e is the energy of the IR photons and the normalization value is fixed at 
£0 = 1' 10~^ eV. One can formally use Eq.(|I]) for an unlimited range of energies 
because at e < 3 • 10^^ eV the microwave radiation dominates while at e > 0.3 eV 
optical radiation does. 

The probability of absorption (the inverse absorption length) of high energy 
photons with energy is given by 0: 

dW 4 foo , , o , n(e) 
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where mg is the electron mass, Ec is the photon energy in the centre of momentum 
system of two colliding photons and cr(ec) is the cross-section for pair production 

-'^The high flux of IR radiation predicted theoretically has recently been questioned by Whipple 
observations m of Mkn421 at energies of a few TeV. 



+ ^ — i> e"*" + e ), which is given in terms of velocity v = {1 — ml/e'lY^'^ as [rU|: 
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where ctt is the Thompson cross-section. 
After simple calculations one obtains: 
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For integer values of u this integral can be solved, e.g., 

$1 = 14/9, $2 = 22/45, $3 = 56/225. 



(4) 
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Over the whole range 1 < z/ < 3, a parametrization accurate to better than one 
percent is given by 



$^ ^ 0.0791 + 1.857e" 
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+ 11.43e 



-2.8761/ 
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For convenience, we tabulate for some selected values of u in Table 1. 
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Table 1: The value of the function of Eq. (|]) for some values of u. 

Now one can put an upper limit on the density of IR photons for the energy 
range 10 — 50 TeV (the range of HEGRA detectability for the set of detectors used 
in the analysis in [|l[]) imposing the condition labs > ^s, where = {c/Hq)z is a 
distance to a HEGRA source^ (we use the Hubble constant Hq = 75 kms~^Mpc~^). 

To explain the HEGRA results we have to assume labs ~ at ^ 50 TeV. 
This condition fixes the flux of IR radiation for a given u. The calculated fluxes are 
exposed in Fig.|l| for different values of u. One can see that these fluxes are much 
lower than ones estimated before. 

The probability of photon absorption (inverse absorption length) is displayed 
in Fig. 2 together with inverse distance to the source 0116+39 (redshift z=0.059) 
shown by the horizontal line. One can observe the dramatic increase of absorption 
at Ej = 100 TeV. The differential flux is suppressed by a factor of 60 — 70 in 
comparison with that at 50 TeV, and integral spectrum is suppressed even more 
due to the sharp cutoff of spectrum at 100 TeV. This may explain why the source 
is not observed by the CASA detector. 



^Since all three HEGRA sources have redshifts z < 0.06, the effects of cosmological evolution 
are very small. 



Photons of the intrinsic radiation with energy between 50 and 100 TeV are ab- 
sorbed on IR radiation and photons with higher energies on microwave radiation. 
The produced electrons and positrons emit photons by scattering off the microwave 
radiation. However, on a scattering length, they are strongly deflected in intergalac- 
tic magnetic field and thus the cascade radiation is spread over a large solid angle. 
This makes the contribution of the cascade radiation to the flux within the angle 



of observation negligible For the energy range of interest, > 1 ■ 10 GeV, 



the energy of radiating electrons is Ee > 2 ■ 10'^ GeV, and on scattering length, 
/ = l/aTUbb, the electron is deflected in magnetic fleld to the angle larger than 
0.01 rad, unless the magnetic field in intergalactic space is smaller than 2 • 10"^^ G, 
which may appear unrealistic. 

Let us consider now neutrino radiation. If the observed gamma radiation is 
produced by protons then the ratio of intrinsic neutrino to gamma ray flux is 1.1 
for 7 = 2.1 and 1.0 for 7 = 2.3. 

Deep underground, the high-energy neutrino flux is accompanied by an equi- 
librium muon flux, i.e the ratio of these fluxes is depth independent. This ratio is 
determined by the neutrino-nucleon interaction and the muon energy losses. For the 
calculations of muon fluxes displayed in Fig. 3, we used the ratios from Ref.0,|]l3 



For values of 7 around 2.1-2.3 the muon spectrum becomes steeper for muon en- 
ergies higher than ~ 1 TeV. This is due to the fact that for energies below 1 TeV 
both the neutrino-nucleon cross section and muon pathlength grow approximately 
linearly with energy, while above this value they grow very slowly. 

As a consequence, most of the muons crossing the detector have energies about 
1 TeV. The calculated flux corresponds to ~ 10 muons with energies higher than 
1 TeV crossing the area 1 km^ per year. They can be marginally detected by the 
future detectors like AMANDA. It should be realised that the number of visible 
neutrino sources could be much larger than the number of TeV gamma sources 
because of the absence of absorption of neutrinos. 

In conclusion, the observations of HEGRA taken at face value imply that the 
flux of extragalactic IR radiation is much lower than was predicted before. The 
combined results of HEGRA and CASA favour the explanation of a 50 TeV cutoff 
due to absorption on IR radiation. In this case the flux of IR radiation is almost 
precisely flxed. The flux at E ^ 1 — 2 TeV is within detectability power of, say, the 
Whipple telescope. The accompaning neutrino radiation can be marginally detected 
by a future 1 km^ detector like an expanded AMANDA. 
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Figure 1: Energy density of the extragalactic radiation field as a function of photon 
energy (logarithmic scales). The curves shown are the cosmic microwave background 
(CMBR), the lowest curve of the analysis of MacMinn and Primack [4], the upper 
limit from the analysis of Biller et al. [10], and the flux needed for interpretation 
of HEGRA data according to our analysis for two values of the slope parameter v 
of the IR flux, V — 2 and v — 2.5. 




Figure 2: Probability of plioton absorption (inverse absorption length) on microwave 
radiation (CMBR) and IR radiation in case oi u — 2 (solid curve ) and u — 2.5 
(dashed curve). The curve labelled z — 0.059 shows the inverse distance to the 
source 0116+319. 
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Figure 3: Deep underground equibrium muon flux produced by high energy neutri- 
nos from 0116+319 source. The energy losses of muons are taken for water. The 
sohd curve corresponds the neutrino spectrum index 7 = 2.3, and dashed curve - 
to 7 = 2.1 Flux of muons for the neutrino spectrum index 7 = 2.3 is higher at low 
energies than that for 7 = 2.1 because both fluxes are normahzed by the HEGRA 
gamma ray flux at £■ = 50 TeV. 



